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The Natural Killer (NK) cell population is composed of subsets of varying sizes expressing 
different combinations of inhibitory receptors for MHC class I molecules. Genes within the 
NK gene complex, including the inhibitory receptors themselves, seem to be the primary 
intrinsic regulators of inhibitory receptor expression, but the MHC class I background is 
an additional Modulating factor. In this paper, we have performed a parallel study of the 
inhibitory receptor repertoire in inbred mice of the C57BI/6 background and in a cohort 
of 44 humans. Deviations of subset frequencies from the "product rule (PR)," i.e., dif- 
ferences between observed and expected frequencies of NK cells, were used to identify 
MHC-independent and MHC-dependent control of receptor expression frequencies. Some 
deviations from the PR were similar in mice and humans, such as the decreased presence 
of NK cell subset lacking inhibitory receptors. Others were different, including a role for 
NKG2A in determining over- or under-representation of specific subsets in humans but not 
in mice. Thus, while human and murine inhibitory receptor repertoires differed in details, 
there may also be shared principles governing NK cell repertoire formation in these two 
species. 
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INTRODUCTION 

The activity of natural killer (NK) cells is controlled by balanc- 
ing inputs from activating and inhibitory receptors. The most 
important ligands for inhibitory receptors are MHC class I mole- 
cules. Because normal cells express high levels of MHC class I, they 
are most often protected from NK cell killing. In contrast, target 
cells expressing downregulated levels of MHC class I are seen as 
"missing self" and killed (Ljunggren and Karre, 1990). Interactions 
between inhibitory receptors and MHC class I also affect NK cell 
development and tunes the responsiveness of the NK cell in an 
"education" process (Hoglund and Brodin, 2010). NK cells sense 
MHC class I in two ways: by direct recognition of individual MHC 
class la alleles or by indirect recognition of an MHC class lb allele 
complexed with peptides derived from class la molecules. Both 
pathways are represented in mice and in humans, suggesting that 
strong selective forces have acted to maintain their parallel use. 
Inhibitory recognition of MHC class lb molecules is performed 
by the NKG2A receptor in mice and humans, where it recognizes 
Qa-l b and HLA-E, respectively. Mice and humans differ, however, 
in the types of inhibitory receptors they use for direct recognition 



of MHC class la molecules. While mice use lectin-like receptors 
of the Ly49 family, human NK cells use killer immunoglobulin- 
like receptors (KIR). In humans, KIR immunogenetics have been 
coupled to autoimmune diseases (Martin et al., 2002), infections 
(Martin et al, 2002; Khakoo et al, 2004; Pelak et al, 2011), and 
outcomes in transplantation (Ruggeri et al., 2002; Cooley et al, 
2009; Venstrom et al., 2012), suggesting that KIR genes take part 
in immunological fine-tuning. KIR-HLA interactions have also 
been proposed to influence NK cell-mediated remodeling of spi- 
ral arteries in the placenta affecting the outcome of pregnancy 
(Hiby et al., 2004). Altogether, these associations imply KIR genes 
as major regulators of human biology. 

Despite the major structural differences between KIR and Ly49 
proteins, the central properties of the two gene families are remark- 
ably similar. First, both families are polygenic, containing between 
15 and 20 genes each (Makrigiannis and Anderson, 2000; Moesta 
and Parham, 2012) and individuals (or mouse strains) vary with 
respect to which of these genes they carry. Secondly, both families 
show polymorphism in individual genes. Thirdly, the expression 
of the gene product at the cell surface is stochastic, giving rise to 
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highly variegated repertoires (Andersson et al, 2009; Brodin and 
Hoglund, 2008). The extent of receptor polymorphism is most 
studied in humans, and new alleles for individual KIR genes are 
discovered continuously (Jiang et al, 2012). Polymorphisms in 
mice are less well characterized since focus so far has been on a few 
inbred mouse strains: C57B1/6 (B6), Balb/C, NOD, and 129 (Mak- 
rigiannis and Anderson, 2000; Patel et al, 2010). This comparison 
has reveled major differences in terms of gene content, but has also 
shown sequence differences suggesting polymorphisms in single 
genes. From those data, a major diversity of NK cell receptors can 
be inferred also in mice, but until the NK gene complex in wild 
rodents from different locations has been analyzed, the extent of 
genetic variation in the mouse remains unknown. 

The magnitude of sequence polymorphisms in individual KIR 
genes is extensive (Jiang et al, 2012). One would think that this 
polymorphism has developed to match the rapid evolution of indi- 
vidual HLA class I alleles, making room for a large variation in 
individual KIR/HLA-interaction patterns. Surprisingly, however, 
the specificity of KIR/HLA interactions seems to be rather con- 
served at the level of HLA recognition, and only four major HLA 
recognition motifs have been identified, each recognized by one 
or several KIRs (Moretta and Moretta, 2004). The role of KIR 
allelic polymorphism in this limited world of ligands is therefore 
unclear. It may affect interaction strength more than overall lig- 
and specificity (Hilton et al, 2012), and could also influence other 
properties of KIR genes, such as expression patterns, that may be 
unrelated to the HLA binding capacity of the receptors (Yawata 
et al., 2008). KIR and Ly49 receptors bind MHC class la molecules 
in a different manner (Achour et al, 1998; Chen et al., 2009). Data 
from mice nevertheless suggest that Ly49 recognition may also be 
flexible in terms of the number of ligands that are recognized by 
each Ly49 receptor (Hanke et al., 1999; Johansson et al., 2009). 
Because of the limited number of individual Ly49 and MHC class 
I alleles that have been studied, it is not possible to know if Ly49 
recognition follows a similar recognition rule as KIR receptors, 
or if recognition is directed against unique residues in individual 
MHC class I alleles. 

A remarkable feature of KIR and Ly49 receptors is the stochas- 
tic and variegated expression patterns that are seen in individual 
mouse and human NK cells. The genetic regulation of this process 
is only partly understood, but the outcome is a "repertoire" of 
NK cells expressing anywhere from zero to five or six inhibitory 
receptors (Yawata et al., 2008; Andersson et al, 2009; Schonberg 
et al., 201 lb; Brodin et al., 2012). Because NK cell repertoires differ 
between individuals, it has been suggested that KIR gene poly- 
morphisms, gene copy number content, epigenetic factors, and 
environmental cues cooperate in repertoire formation (Uhrberg, 
2005; Pascal et al, 2006). The presence of MHC class I ligands for 
inhibitory receptors also affect the frequencies of individual NK 
cell subsets, which has been interpreted in terms of an "adapta- 
tion" process to ensure optimal "missing self" recognition in each 
individual (Hoglund and Brodin, 2010). The evidence for this 
is strongest in mice, where direct comparisons of the inhibitory 
receptor repertoires in MHC class I-deficient and sufficient mice 
have been made (Held et al., 1996; Salcedo et al., 1998; Fahlen 
et al, 2001; Brodin et al, 2009, 2012). In humans, links between 
the presence of self HLA ligands and expansion or contraction of 



certain NK cell subsets have been difficult to demonstrate (Shilling 
et al, 2002; Andersson et al, 2009). High resolution KIR pheno- 
typing has however shown effects of HLA on the KIR repertoire 
(Yawata et al., 2006; Yu et al., 2007). Our recent data suggest that 
KIR ligands do have a clear effect on the KIR repertoire but only 
in individuals infected by CMV (Beziat et al, 2013). 

In this article, we will discuss the formation of murine and 
human inhibitory receptor repertoires. We will first describe the 
repertoire in B6 mice, and examine how the observed repertoires 
differ from the expected repertories under the assumption that 
repertoire composition is completely random. Using this analysis, 
we identified both MHC-independent and dependent effects of 
the repertoire. We will finally compare the murine repertoire of 
one single strain with that of a cohort of humans homozygous for 
the group A KIR haplotype. 

MATERIALS AND METHODS 
MICE 

Inbred mice of the B6 genetic background were used. B6 mice 
lacking functional MHC class I heavy chain genes (B6.K b D w - 
mice) were a kind gift from Francois Lemonnier, Institut Pasteur, 
France. B6.K b D b_/_ mice were obtained by crossing mice made 
deficient for the K b and D b genes respectively (called B6.K b ~'~ and 
B6.D b ~'~ mice) as described previously (Perarnau et al., 1999; 
Lemonnier, 2002). To generate mice expressing D d as the only 
MHC class la allele, H2D d transgenic B6 mice (Hoglund et al, 
1988) were crossed to B6.K b D b ~ / ~ mice, followed by further back- 
cross of the progeny to B6.K b D b ~ / ~ mice. Individuals homozygous 
for the K b and D b mutations and positive for the D d transgene 
were identified using FACS analysis. A similar mating scheme was 
followed to generate mice expressing L d alone starting with H2L d 
transgenic B6 mice (Johansson et al., 2005, 2009) . All animal exper- 
iments were performed in accordance with relevant institutional 
and national guidelines and regulations and conform to relevant 
regulatory standards. Ethical evaluations were performed by the 
animal ethics committee of northern Stockholm. Mice in each 
experiment were age (6-10 weeks) and sex matched. Both female 
and male mice were used. 

HUMAN SUBJECTS 

Experiments on human blood samples were approved by the 
regional ethics committee (Stockholm, Sweden, approval num- 
ber 2006/229-31/3). Buffy coats were prepared from peripheral 
blood of healthy human donors (ages approximately 20-60 years 
old) and separated by density gradient centrifugation (Ficoll- 
Hypaque™, GE Healthcare Bio-Sciences AB, Uppsala, Sweden). 
The cohort of humans used in this study is the same as the one used 
in our previous study, and includes 44 individuals homozygous for 
the group A KIR haplotype (Andersson et al, 2009). 

FLOW CYTOMETRY 

Murine NK cells were isolated from the spleen and human NK 
cells from peripheral blood (Andersson et al, 2009; Brodin et al, 
2012). The antibody panels and protocols for multi-color flow 
cytometry used for the repertoire studies in mice and humans 
have been described in detail (Fauriat et al., 2008; Anders- 
son et al, 2009; Brodin et al, 2012). In brief, human NK cell 



Frontiers in Immunology | NK Cell Biology 



March 2013 | Volume 4 | Article 65 | 2 



Sternberg-Simon et al 



Human and murine inhibitory receptor repertoires 



receptor antibodies were KIR3DL1 (DX9), KIR2DL1/S1 (EB6), 
KIR2DL2/3/S2 (GL183), KIR3DL2 (DX31), and NKG2A (zl99). 
For the mouse analysis, Ly49A (YE1/48), Ly49C (4L0331 1), Ly49I 
(YLI-90), Ly49G2 (4D11), and NKG2A/C/E (20d5) antibodies 
were used. Gating strategies for the FACS data used here have been 
described (Fauriat et al., 2008; Brodin et al, 2012). In summary, 
murine NK cells were identified as NK1.1 + CD19— CD3— and 
human NK cells as CD56 + CD14-CD3-. All stainings in both 
species excluded dead cells and doublets. Receptor expression was 
evaluated on these gated subsets. 

PRODUCT RULE CALCULATIONS 

The observed exclusive expression probabilities of combinations 
of zero to five receptors were calculated from the experimental 
data. In addition, we calculated the expected exclusive expres- 
sion probabilities for these combinations of receptors, based on 
the overall expression frequency of each receptor in the data, and 
assuming that the expression follows the product rule (PR). Devi- 
ations of the observed from the expected frequencies of expression 
of receptor combinations are plotted as log2 (observed/expected). 
In box plot figures, the box signifies upper and lower quartiles, and 
the median is represented by a black line within the box. Bottom 
and top whiskers represent the 5 and 95% values, respectively. 

RESULTS 

THE INHIBITORY RECEPTOR REPERTOIRE OF B6 MICE 

A flow cytometry panel was used to simultaneously detect the 
frequency and expression levels of the following five inhibitory 
receptors on NK cells in B6 mice: Ly49A, Ly49C, Ly49G2, Ly49I, 
and NKG2A. Using this antibody panel, we determined the fre- 
quencies of NK cells expressing all possible combinations of these 
receptors, in total 32 subsets. Six experiments were performed, 
including 37 individual mice. All 32 possible subsets were present 



in the repertoire but in greatly varying frequencies (Figure 1). Of 
note, NK cells lacking all five receptors were among the most preva- 
lent, reaching 7.4% in this analysis (Figure 1). The standard devia- 
tions were rather large for most subsets, which was surprising given 
the genetic homogeneity of this cohort of inbred mice (Figure 1). 
This variation was to a large extent due to inter-experimental vari- 
ation of subset frequencies (Figure Al in Appendix). However, 
there was also intra-experimental variation, of similar size in each 
experiment, that could reflect a biological variation, transient or 
permanent, in the expression of these inhibitory receptors. 

OBSERVED SUBSET FREQUENCIES DEVIATE FROM THOSE EXPECTED 
FROM PRODUCT RULE CALCULATIONS 

In a model in which expression of individual genes is indepen- 
dently regulated, the fractions of cells expressing combinations of 
genes can be predicted using "PR" calculations. According to the 
PR, the joint probability of expression of two or more receptors 
is given by the product of the individual expression frequencies 
of these receptors (Mehr et al, 2012). Calculations can be per- 
formed in two ways: "inclusive" calculations examine the joint 
expression probability of a combination of receptors, regardless 
of other receptors expressed by the same cell. In contrast, "exclu- 
sive" calculations examine the joint expression probability of a 
combination of receptors, excluding cells that express any other 
receptors (Figure 2A). In this study, we have used the "exclusive" 
way of calculating the PR, since it predicts not only frequencies of 
NK cells expressing multiple receptors, but also subsets without 
receptors or with single receptors (Mehr et al., 2012). By com- 
paring the observed and predicted expression frequency of each 
receptor combination, inter-dependencies between receptors and 
specific deviation in different genetic backgrounds, such as the 
comparison between MHC-deficient and B6 mice, can be revealed. 
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FIGURE 1 | The NK cell inhibitory receptor repertoire in B6 mice. The 

observed frequency of each of 32 possible combinations, averaged for 37 
mice, analyzed in six independent experiments. Error bars indicate the 



standard deviations. The lower panel indicates, for each combination, whether 
an individual receptor is expressed in this combination or not (black and white 
squares, respectively). 
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FIGURE 2 | The observed repertoire in B6 mice deviates from the product 
rule. (A) The inclusive (left) and exclusive (right) calculation of the probability to 
express two out of three possible receptors. If the probability to express each 
of the receptors A, B, and C is P(A), P(B), and P(C), respectively, then the 
inclusive probability to express A and B is P(AB) = P(A) x P(B) and it includes 
all the cells that expresses both A and B, regardless to the expression of C, 
and therefore includes combinations AB and ABC. The exclusive probability to 



express A and B is P(AB) = P(A) x P(B) x [1 -P(C)], and it includes only cells 
that express A and B, but do not express C. (B) The observed (black bars) and 
expected (white bars) repertoire of one representative mouse. (C) Deviations 
from the product rule in 37 B6 mice, shown as box plots of log 2 (observed 
frequency/expected frequency). Negative values indicate that the respective 
combination is under-expressed, compared to the expected frequency, and 
positive values indicate over-expression. 



The observed and expected frequencies of all 32 Ly49 subsets 
in B6 mice are exemplified by a representative mouse (Figure 2B), 
as well as a summary of all mice in the cohort plotted as Log2 
ratios between observed and expected frequencies (Figure 2C). 
Deviations were common and occurred in both directions, i.e., 
frequencies were both increased and decreased relative to the 
expected (Figure 2C). Examination of subsets expressing single 
Ly49 receptors for self MHC class I in B6 mice, Ly49C, Ly49I, or 



NKG2A (Brodin et al., 2012), did not reveal a consistent devi- 
ation pattern. In fact, these three subsets behaved completely 
differently; the Ly49C subset was over-represented, the Ly49I 
subset under-represented, and the NKG2A subset was close to 
unchanged. Subsets expressing two or three receptors showed 
deviations in both directions while subsets expressing no, four, 
or five receptors were in most cases less frequent than expected 
(Figure 2C). 
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FIGURE 3 | Deviations from the product rule of single receptor 
combinations in B6 and MHC mice. (A) The observed frequencies of NK 
cells expressing the indicated single receptor combination in 43 MHC and 
37 B6 mice (left and right panel, respectively) are plotted against the 
expected frequency, under the product rule assumption. A perfect fit to the 
product rule is shown as a straight line. (B) Deviations from the product rule 
in 37 B6 mice, after subtracting, for each mouse and for each receptor 
combination, the average deviation measured in 43 MHC mice for the same 
combination. 



PRODUCT RULE DEVIATIONS ARE CAUSED BY A COMBINATION OF 
MHC CLASS l-DEPENDENT AND INDEPENDENT EVENTS 

To investigate to what extent PR deviations of the three sub- 
sets expressing single self-specific receptors were MHC class I- 
dependent, we performed the same analysis in a cohort of 43 mice 
lacking MHC class I molecules. Interestingly, deviations of these 
subsets were seen also in MHC-deficient mice (Figure 3A, left 
panel). For NK cells expressing only Ly49C, most MHC-deficient 
mice contained fewer such NK cells than predicted by the PR, but 
some mice showed an over- representation (Figure 3A, top left). 
This pattern was different from that in B6 mice, where all mice con- 
tained a higher frequency of these cells than expected (Figure 3A, 
top right). The pattern for Ly49I (under- representation) was very 
similar in MHC-deficient and in B6 mice (Figure 3, middle), 
while NKG2A single-positive NK cells were more frequent than 
expected in all MHC-deficient mice, yet showed a close to expected 
expression pattern in B6 mice (Figure 3A, bottom). 

To reveal the total pattern of MHC class I-driven deviations 
from the PR of all 32 subsets in B6 mice, we subtracted the aver- 
age deviation of each subset in the cohort of MHC-deficient mice 
(shown in Figure A2 in Appendix) from the deviation in each 
individual B6 repertoire. The resulting Log2 plot (Figure 3B) 
showed several differences compared to the unsubtracted B6 data 
(Figure 2C). A general effect was an increase in the number of 
subsets which were over- represented compared to expected, a pat- 
tern that was particularly clear for subsets expressing two receptors 
(Figure 3B) . However, there were also subsets that went from being 
slightly over-represented to under-represented, such as NK cells 
expressing only NKG2A. 

DISTINCT REPERTOIRE PERTURBATIONS IN MICE EXPRESSING SINGLE 
MHC CLASS I ALLELES 

To test whether different MHC class I alleles produced distinct 
changes of the repertoire, we compared four strains expressing 
single MHC class I alleles (Johansson et al, 2005, 2009) and sub- 
tracted, from their PR deviations, the average PR deviations of 
each receptor combination in MHC-deficient mice (Figure 4). 
NK cells lacking all inhibitory receptors were fewer than expected 
in all strains, suggesting that this is a common deviation induced 
by all MHC class I alleles. Secondly, H2D d promoted a devia- 
tion of the Ly49A single receptor-expressing subsets to an extent 
not observed by the other alleles (Figure 4), which is in line 
with our previously published results (Brodin et al., 2012). On 
the same note, Ly49C single-positive NK cells were more preva- 
lent than expected in H2K b mice compared to H2D b mice, in 
line with the identification of H2K b as a strong Ly49C ligand in 
B6 mice (Hanke et al, 1999). Thirdly, H2D b and H2L d , which 
are very similar in terms of sequences and expression proper- 
ties, produced a nearly identical perturbation pattern with several 
deviations, including a distinct over-representation of one par- 
ticular subset expressing the combination of Ly49A, Ly49I, and 
NKG2A (Figure 4, middle two panels). Finally, we noted that the 
MHC class I allele with the broadest specificity of the ones tested, 
H2D d (Hanke et al, 1999; Johansson et al, 2005, 2009) produced 
a deviation pattern that was most similar to the one seen in B6 
mice, where two alleles exert an educating impact on the NK cell 
repertoire. 



THE AVERAGE HUMAN REPERTOIRE IS HIGHLY DIVERSE 

We next performed a similar repertoire analysis in a cohort of 
44 healthy donors that were homozygous for the group A KIR 
haplotype (Fauriat et al, 2008; Andersson et al, 2009). Superfi- 
cially, the average human KIR repertoire was quite similar to the 
mouse Ly49 repertoire, with more than 1 0% of NK cells expressing 
no inhibitory receptor, and an overall larger fraction of subsets 
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FIGURE 4 | MHC-dependent deviations from the product rule in single 
MHC mice. The deviations from the product rule of 37 K b mice, 6 D b mice, 
6 L d mice, and 34 D d mice, after subtracting, for each mouse and each 
receptor combination, the average deviation observed in 43 MHC mice for 
the same combination. 



expressing one or two inhibitory receptors compared to subsets 
expressing more than three receptors (Figure 5A). As expected 
given the genetic heterogeneity of outbred humans, we observed 
high variation between individuals, resulting in large standard 
deviations. The nature of the variation became clear when four 
individual human repertoires, expressing an identical setup of 
KIR ligands (not shown), were plotted side-by-side (Figure 5B). 
Among these four donors, the first (donor #5) had an inhibitory 
receptor repertoire that resembled that of B6 mice in terms of the 
distribution of subsets expressing distinct receptor combinations. 
The other three (donors # 19, 25, 43) were different, and showed 
either a more restricted expression pattern (donors #19 and 43) or 
a more equal distribution of subsets (donor #25). Altogether, these 
data support previous results showing substantial KIR ligand- 
independent variations in the inhibitory receptor repertoires in 
humans (Gumperz et al., 1996; Shilling et al, 2002; Andersson 
etal.,2009). 
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FIGURE 5 | The human repertoire is highly diverse. (A) The observed 
frequency of each of 32 possible combinations, averaged for 44 haplotype 
A donors. Error bars indicate the standard deviation. (B)The observed 
repertoires in four individual haplotype A donors. 



PATTERNS OF PRODUCT RULE DEVIATIONS REVEAL DIFFERENCES 
BETWEEN MOUSE AND HUMAN REPERTOIRES 

Because PR calculations are performed similarly in mice and 
humans, we analyzed the repertories in these species side-by-side 
by grouping subsets according to the number of receptors they 
express. Observed and expected frequencies of each subset were 
then visualized in a dot plot. Four observations were made in this 
comparison. First, NK cells devoid of inhibitory receptors were less 
frequent than expected in both B6 mice and in outbred humans 
(Figure 6, upper two plots), suggesting that under- representation 
of such NK cells is a general rule. Secondly, for NK cells express- 
ing single inhibitory receptors, NKG2A-expressing NK cells were 
strongly over-represented in humans (Valiante et al, 1997; Fau- 
riat et al, 2008; Yawata et al, 2008) while no deviation of this 
subset was seen in B6 mice. Linked to this notion, for sub- 
sets expressing three to five receptors, the presence or absence 
of NKG2A determined whether or not the subset was over- or 
under- represented relative to the expected frequency (Figure 6). 
A similar role for the NKG2A receptor was not observed in B6 
mice. Thirdly, in humans, three donors showed a strong over- 
representation of KIR2DL2/3 -positive NK cells (red dots) and one 
donor showed a similar over-representation of KIR2DL1 -positive 
NK cells (black dots). It is tempting to speculate that these devi- 
ations resulted from a CMV-driven and KIR ligand-dependent 
skewing of the repertoire, as we recently noted in a much larger 
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FIGURE 6 | Deviations from the product rule in B6 mice and in humans. panel) and 44 A haplotype human donors (right panel). A perfect fit to the 
The observed frequency of NK cells co-expressing zero to five receptor product rule is shown as a straight line. Each symbol indicates a specific 

combinations are plotted against the expected frequency, in 37 B6 mice (left receptor combination, as shown on the panels beside the figures. 



cohort of donors (Beziat et al., 2013). However, this remains a 
speculation since serological data regarding previous exposure to 
CMV is lacking in this cohort. Finally, subsets expressing all five 
receptors were strongly under-represented in B6 mice but not 
in humans. Thus, our comparison, despite suffering from some 
drawbacks (see Discussion) revealed similarities as well as differ- 
ences between B6 mice and humans in the formation of inhibitory 
receptor repertoires. 



DISCUSSION 

A distribution of cellular subsets based on variegated expression of 
inhibitory receptors is a unique and intriguing property of human 
and murine NK cells (Moesta and Parham, 2012). To understand 
why this expression pattern has evolved, it is important to con- 
sider that the number of different inhibitory receptors an NK 
cell expresses determine its function under certain circumstances 
(Hanke and Raulet, 2001; Yu et al., 2007). For example, to elicit 
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a "missing self" response, an NK cell expressing two different 
inhibitory KIR for self MHC class I will require downregulation of 
the HLA ligands for both of those. In contrast, an NK cell carrying 
a single Ly49 or KIR receptor will respond to targets missing single 
HLA class I alleles and will thus be sensitive to specific alterations 
in HLA expression. 

A stochastic and independent expression pattern of each KIR 
or Ly49 gene represent one solution to generate NK cells with 
restricted numbers of Ly49 or KIR receptors (Cichocki et al, 201 1; 
Gays et al, 2011). In this process, the development of NK cells 
expressing several or zero receptors is unavoidable, however, lead- 
ing to the formation of subsets that might either be constantly 
inhibited by self MHC (many self receptors), or potentially autore- 
active as they are not inhibited by self (no receptors). Subsets 
lacking inhibitory receptors for self may cause harm, but there 
are mechanisms to make them tolerant (Fernandez et al., 2005; 
Johansson et al, 2005; Kim et al, 2005; Anfossi et al., 2006). In this 
general respect, NK cell education bears similarities to T cell selec- 
tion, in which harmful T cells may arise that are kept in control 
by complementary mechanisms. Which NK cell functions in vivo 
that are controlled by the "missing self" principle is unknown (Orr 
et al, 2010) and this question represents one area of intensive 
research. 

Because of the genetic homogeneity and standardized hous- 
ing conditions of inbred mice, genetic factors that have been 
implicated in the control of KIR expression in humans (such 
as promoter use, gene copy number, allelic variants, and infec- 
tion history) are neutralized. It was therefore surprising that the 
observed repertoires showed a rather large variation between indi- 
vidual mice. This variation was caused by a combination of inter- 
experimental variation (which was largest) andintra-experimental 
variation. It is not known if this is caused by technical aspects of 
the analysis, or if individual mice differ from each other due to bio- 
logical fluctuations in receptor expression. Longitudinal analysis 
of NK cell repertoires by repeated analyses of individual mice may 
give some clues to this question. 

The significant and consistent deviations from the PR in mice 
deficient for MHC class I heavy chains was unexpected. One expla- 
nation for these deviations is that the expression probability of the 
various genes in the NKC are not independent, and that the expres- 
sion of one gene in NKC may affect the probability of another gene 
to be expressed. Such dependencies in the acquisition of multi- 
ple inhibitory KIRs has been suggested in humans to explain the 
extensive PR deviations observed in NK cells expressing three and 
four inhibitory KIRs (Gardiner et al., 2001; Andersson et al, 2009). 
Local access to transcription factors or simultaneous demethyla- 
tion of multiple promoter regions, could also play a role. It is also 
possible that ligands for inhibitory receptors other than classical 
MHC class I molecules may affect receptor skewing. Such ligands 
can at least influence licensing of subsets, as shown recently for 
the non-classical H2-M3 gene (Andrews et al, 2012). Another 
possibility is that there are biochemical dependencies between 
inhibitory receptors that affect their transport or coexpression at 
the cell surface. 

A remaining question, especially in humans, is whether the 
repertoire is determined only by genetically controlled mecha- 
nisms acting intrinsically within the NK cell, or if there are aspects 



of the repertoire that are subject to selection to optimize it for 
missing self recognition. For example, are there processes that 
positively select certain subsets and negatively select other subsets, 
depending on the MHC class I setup? When testing this ques- 
tion, MHC class I-independent skewing effects must be taken into 
account. In our study, we decided to subtract the deviations in 
MHC-deficient mice to allow a better resolution of the MHC class 
I-controlled repertoire in B6 mice. By doing so, we were able to 
draw some general conclusions, as follows. First, NK cell subsets 
that express no inhibitory receptors (out of those studied) were 
under-represented in B6 mice. These are potentially harmful as 
they may not be inhibited by MHC class I-expressing self cells. 
Their under-representation adds an additional layer of protection 
to the hypo-responsiveness that such cells exhibit in "missing self" 
contexts (Fernandez et al, 2005; Kim et al., 2005; Anfossi et al, 

2006) . 

Secondly, a large fraction of subsets were more frequent 
than expected in B6 mice, including all subsets expressing two 
inhibitory receptors, and most (but not all) subsets expressing 
three receptors. This finding suggests that MHC class I-mediated 
effects are broader than previously thought. Thirdly, NK cells 
expressing one known self receptor in isolation, Ly49C, Ly49G2, 
and NKG2A, showed different skewing patterns. This latter finding 
speaks against our previously proposed notion of positive selection 
as a general mechanism to enrich for NK cell subsets expressing 
a single self receptor. Still, data on the selective accumulation of 
Ly49C in the presence of H2K b (this paper) and Ly49A in the 
presence of H2D d (this paper and Brodin et al., 2012) suggest that 
positive selection may operate for some such subsets, even if other 
subsets with expression of a single self receptor behave differently. 
One interesting property of Ly49A and Ly49C, but not Ly49I, that 
comes to mind in this respect is their abilities to bind their respec- 
tive MHC class I ligands in cis (Kase et al., 1998; Doucey et al, 
2004; Andersson et al., 2007; Back et al., 2007; Scarpellino et al, 

2007) . NKG2A has not been studied in this respect. Finally, subsets 
expressing four or five receptors tended to be under-represented in 
B6 mice. Overall, the difference in repertoire composition between 
B6 mice and MHC-mice confirm a role for self MHC class I 
in shaping the repertoire. Much more work is required before 
we understand all alterations, in particular alterations of subsets 
expressing two or three receptors. At this point, we conclude that 
the effects of MHC class I was surprisingly broad and involved 
most subsets. 

In humans, it has been more difficult to detect HLA class I- 
mediated skewing of the KIR repertoires. In fact, because random 
acquisition of KIRs in itself results in repertoires dominated by 
NK cells expressing few KIRs (most of which express at least one 
self-specific KIR), the existence of an intrinsic selection process 
has been questioned (Andersson et al, 2009). Nevertheless, some 
studies have detected increased frequencies of NK cells express- 
ing receptors for HLA-C (Schonberg et al., 201 lb) and HLA-Bw4 
(Yawata et al., 2006) in individuals expressing ligands for these 
receptors. In our previous work, we failed to see a global skewing 
toward self (Andersson et al., 2009) that was predicted by models of 
selection favoring NK cells with self-specific receptors. Although it 
is possible that the inherent variation in KIR expression between 
individuals is too large to detect global effects, the finding that 
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neonatal NK cells lack even subtle skewing toward self (Schonberg 
et al, 2011a), suggest that human repertoires are formed ran- 
domly and independently of HLA class I. The dichotomy between 
subtle skewing effects in human adults but nor in cord blood 
may be explained by our identification of CMV as a cofactor in 
HLA-driven repertoire skewings (Beziat et al, 2013). 

The mechanisms underlying MHC class I-associated pertur- 
bations of the repertoire are unknown. Because all subsets in the 
repertoire add up to 100%, changes can be explained either by pro- 
liferation of some or death of other subsets. Alternatively, MHC 
class I may impose effects on the probability to express an addi- 
tional receptor to the ones a developing cell already has, without 
effect on death or cell division. Our previous data showed that 
Ly49A + NK cells from mice expressing the Ly49A ligand H2D d 
responded better to IL-15 than Ly49A+ NK cells from MHC- 
deficient mice (Brodin et al, 2012), suggesting that skewing may 
be mediated by specific expansion of educated NK cells. As already 
mentioned, the fact that individual single MHC class I alleles pro- 
duce specific changes in the repertoire consistent with their known 
specificities for Ly49 receptors, imply positive selection of at least 
some useful subsets as one driver for repertoire adaptations. Recent 
data from humans suggest a similar scenario, at least in the context 
of viral infection. Thus, acute and latent CMV infection is associ- 
ated with a profound imprint in the human KIR repertoire caused 
by an expansion of NK cells expressing self-specific KIRs (Delia 
Chiesa et al., 2012; Foley et al, 20 12; Beziat et al, 20 13). In vitro, NK 
cells educated by self-specific inhibitory KIRs expand specifically 
when stimulated by target cells expressing HLA-E or fibroblasts 
infected with CMV (Beziat et al, 2012, 2013; Charoudeh et al, 
2013). Although it is still unclear whether the expansion of NK 
cells expressing self-specific KIRs is useful for the host or not, the 
resolution of CMV reactivation in a patient with T and B cell 
SCID correlated closely with the expansion and contraction of 
NKG2C+ NK cells expressing self-specific KIRs (Kuijpers et al, 
2008). 

Product rule calculations represent a tool to estimate devia- 
tions of subsets frequencies relative to the expected, assuming that 
individual receptors are independently expressed. To make this 
comparison, the algorithm must be fed the individual expression 
probabilities of all receptors under study. Because these probabili- 
ties are unknown, it is instead assumed from the overall frequency 
of a given receptor in the repertoire. While this procedure makes 
PR calculations possible in all individuals, it also contains a circular 
loop in that the overall frequency of a given receptor will be affected 
by deviations of other subsets, which is what PR calculations intend 
to measure. This caveat carries a risk of underestimating small sub- 
set deviations. On the other hand, deviations observed despite this 
may perhaps be more safely interpreted, even if this must still be 
done with caution, and with the possibility of other potentially 
contributing mechanisms in mind. 

Exclusive calculations of product rule deviations were used to 
analyze an inbred mouse strain and a cohort of humans in parallel. 
Of note, our mouse cohort represent 37 independent observations 
of one single genotype while the human cohort included 44 out- 
bred donors. General conclusions regarding the diversity of the 
repertoire in the population may therefore be made in humans, 
but not in mice. In contrast, the genetic heterogeneity of humans 



may hide specific effects in individuals, while the mouse analysis is 
stronger on that point. To conceptualize the difference, one could 
view the mouse data as an example of a repeated analysis of a single 
case in the human cohort. 

A side-by-side view still gave several insights. Taking a look 
at the human cohort, we first concluded that NK cells lacking 
inhibitory receptors were under-represented. This observation is 
consistent with a model in which MHC class I proteins posi- 
tively select cells expressing inhibitory receptors. Our data on 
B6 mice as well as four single MHC class I mice support this 
as a mechanism also in mice. The most notable feature of the 
human repertoire was the effects of NKG2A expression on NK 
cells, which was manifested as an over-representation of single- 
NKG2A-expressing NK cells and an under-representation of sub- 
sets expressing NKG2A together with KIR. One way to interpret 
this result is to postulate that the expression of NKG2A would 
make subsequent KIR expression less likely. In this scenario, an 
NK cell that express NKG2A as its first inhibitory receptor would 
tend to be over-represented and NK cells expressing NKG2A in 
association with KIR would be less frequent. Perhaps this model 
could also explain why NKG2A has been observed to "buffer" the 
repertoire in some donors (Yawata et al, 2008; Andersson et al, 
2009). Thus, in cases where the probability to express NKG2A is 
high, most NK cells would lack KIR and express only NKG2A, 
while donors with a low expression probability of NKG2A would 
have a repertoire dominated by KIR expression (Valiante et al, 
1997). Intriguingly, in B6 mice, a similar role for NKG2A was not 
seen in the repertoire. As alluded to, it is not known if this is a 
feature of the B6 background or a general property of all mouse 
strains. 

The inverse correlation between NKG2A and KIR in the NK 
cell repertoire may also reflect population of cells at discrete stages 
of differentiation (Bjorkstrom et al., 2010; Beziat et al., 2011). 
NKG2A is expressed by all CD56 bri s ht cells and is gradually lost 
during terminal NK cell differentiation within the CD56 dlm stage, 
correlating with the acquisition of KIRs and expression of CD57 as 
well as loss of CD62L and NKp30 (Bjorkstrom et al, 2010; Juelke 
et al., 2010). The time frame for these events remains unknown. It 
is possible that cells acquire their KIR repertoire very early during 
development and that the original set up of expressed KIRs deter- 
mine whether the cells are "allowed" to lose NKG2A during their 
continued differentiation. 

The term "NK cell education" has been defined by many to 
mean the acquisition of functional competence in NK cells, result- 
ing from interactions with MHC class I molecules. We have pro- 
posed that alterations in the inhibitory receptor repertoire may be 
included in a broader definition of NK cell education (Hoglund 
and Brodin, 2010). Our data, from this and other papers, on reper- 
toire perturbations in single MHC class I mice and CMV-positive 
humans, support this view. The link between functional education 
and repertoire skewing, both in terms of regulation of global NK 
cell responses and in terms of signaling pathways involved remains 
to be determined. In addition, the process of NK cell education is 
still poorly defined in space and in time, and may involve a combi- 
nation of developmental influences on immature precursors and 
tuning of mature NK cells in the periphery. Presently, we favor 
the view that MHC class I molecules educate the NK cell system 
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by many means, one being to license them to become functional 
and the other being shaping the repertoire toward the expression 
of certain combinations of receptors. Further work in our labo- 
ratories is directed toward the understanding of how these two 
processes are linked. 
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APPENDIX 



Variability between experiments with B6 mice 
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FIGURE Al | Repertoire variations in each experiment separately 
(#1-#6) and in all mice together (All Exp). The average frequencies of 
each of the 32 NK cell subsets were first calculated in each group of mice. 
The coefficient of variance (CV, standard deviation/average) for each subset 
was determined and the average of all subsets in each group was plotted 
with the standard deviation. The number of mice in each experiment is 
indicated in parentheses. The CV's were (from left to right): 015, 0.11 , 0.109, 
0.181, 0.182, 0.2, and 0.33. The higher the CV, the higher the variation is in 
the group. 
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FIGURE A2 | Deviations from the product rule in 43 MHC-deficient mice, shown as box plots of log 2 (observed frequency/expected frequency). 

Negative values indicate that the respective combination is under-expressed, compared to the expected frequency, and positive values indicate over-expression. 
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